/ / ' 

NASA Contractor Report 201 739 
ICASE Report No. 97-45 



Divergence Boundary Conditions for Vector Helmholtz 
Equations with Divergence Constraints 


Urve Kangro 
Roy Nicolaides 


NASA Contract No. NASI -19480 
September 1997 

Institute for Computer Applications in Science and Engineering 
NASA Langley Research Center 
Hampton, VA 23681-0001 

Operated by Universities Space Research Association 



National Aeronautics and 
Space Administration 


Langley Research Center 

Hampton, Virginia 23681-0001 



Divergence Boundary Conditions for Vector 
Helmholtz Equations with Divergence Constraints 


Urve Kangro* 

Department of Mathematical Sciences 
Carnegie Mellon University 
Pittsburgh, PA 15213 


Roy Nicolaidest 

Department of Mathematical Sciences 
Carnegie Mellon University 
Pittsburgh, PA 5213 


Abstract 

The idea of replacing a divergence constraint by a divergence boundary 
condition is investigated. The connections between the formulations are 
considered in detail. It is shown that the most common methods of using 
divergence boundary conditions do not always work properly. Necessary 
and sufficient conditions for the equivalence of the formulations are given. 


•This work was supported by the United States Air Force Office of Scientific Research 
under grant F-49620- 94- 1-0311. 

^This work was supported by the United States Air Force Office of Scientific Research 
under grant F-49620- 94- 1-0311. This work was also supported in part by the National Aero- 
nautics and Space Administration under NASA Contract NASI- 19480 while the author was 
in residence at the Institute for Computer Applications in Science and Engineering (ICASE). 


l 




1 Introduction 


A well known problem in computational electromagnetics is the appearance of 
“spurious modes” arising from incorrectly imposing divergence constraints. To 
avoid imposing these constraints directly, which is difficult, it has been suggested 
that they need to be enforced only on the boundary of the domain ([5], [7]). This 
has some significant advantages. For instance, the vector Helmholtz equation 
with a divergence constraint on the field can then be solved using standard finite 
element spaces instead of more complex spaces of edge elements. Furthermore, 
as we will show below, the divergence boundary condition can be treated as a 
natural boundary condition. 

It is essential to know when the interior and boundary formulations of diver- 
gence constraints give the same result and that is the main aim of this work. We 
shall show that the formulations are equivalent if and only if the scalar Poisson 
equation with any smooth right hand side and Dirichlet boundary conditions 
has a solution in H 2 . When this is not true, a simple minded application of 
divergence boundary conditions as in [5] and [7] will normally give incorrect 
solutions with a nonzero divergence. 

Our results for systems with the interior divergence constraint are obtained 
by formulating the problem in the space V defined in the next section. This 
approach avoids enforcing the divergence constraint explicitly — - it is automati- 
cally satisfied. We provide a proof of this fact for interior problems; for exterior 
domains see [4]. Unfortunately, H 1 is not dense in V in general. To use stan- 
dard H 1 finite elements we impose the divergence condition on the boundary 
and formulate the problem in the space C H 1 defined below. We discuss 
the relationship between the V and H^ t formulations in some detail. 

To avoid relatively unimportant technical issues we will give proofs of these 
results for the vector Poisson equation subject to tangential boundary conditions 
and a divergence constraint. In the end we will show how they may be extended 
to other problems including the vector Helmholtz case of electromagnetics. 

The next section gives the strong forms of the governing equations and the 
reasoning behind the use of divergence boundary conditions. It also contains a 
formula for the boundary divergence, which is used to reformulate the divergence 
boundary condition as a natural boundary condition in H 1 setting. Section 3 
sets up a weak form for the equation with the interior divergence constraint 
and proves the coercivity for the weak form, which follows from a compact 
embedding result for vector fields. Sections 4 and 5 prove coercivity for the weak 
form which uses divergence boundary conditions. This time, coercivity follows 
from a close and somewhat surprising connection with the previous weak form. 
This result is based on a trace theorem proved in section 5. 

The coercivity results of sections 3 5 are used in section 6 to prove equiv- 
alence of the strong and weak formulations (Theorem 2). We also discuss in 
what sense the boundary divergence condition is satisfied and prove a formula 
for calculating the divergence on the boundary (Lemma 10 in section 6). Section 
7 contains a discussion of when one can say a priori that the interior and bound- 
ary divergence formulations are equivalent, and when they are not equivalent 
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(i.e. may have different solutions for the same data). It also contains an exam- 
ple showing nonequivalence of the two formulations for nonconvex polyhedra. 
In section 8 we give the modifications to include the Helmholtz equation and 
briefly mention some other extensions. Section 9 mentions some implications 
for numerical approximations of the problem. 


2 Formulation of the problems 

Assume that fi is an open bounded subset of R. V , N = 2 or 3, with a connected, 
Lipchitz boundary I\ 

We will consider the relationship between the problems 


and 



We assume that f E L 2 (Fl) N and the compatibility condition divf = 0 is 
satisfied. Superficially it seems that the two problems are equivalent. Indeed, 
if u is a solution of (1), it is clearly a solution of (2). Conversely, if u is a 
solution of (2), then taking the divergence of the first equation of (2) we get 
that A(divu) = 0, and since divu = 0 on T, one might argue that divu = 0 
in fl, so u is a solution of (1). This is basically the reasoning in [7]. But it 
might not always be true. The difficulty is caused by insufficient regularity 
of the solution. Laplace’s equation with Dirichlet boundary conditions has a 
unique solution in /^(Q), but in general we cannot assume that div u e 
I t may be only in L 2 (Fl) and the homogeneous equation for divu may have a 
nonzero solution in L 2 (fJ) (this happens, for example, in the case of domains 
with reentrant corners). 

In order to discuss the the relationship between (1) and (2) we first have to 
make explicit in what function spaces we look for the solutions. It turns out 
that the natural space for (1) is 

V = {u E L 2 (Q) n | divu E L 2 (tt), curlu E L 2 (Fl) N \ u x n|r = 0}, 

where TV 7 = 3 if N = 3 and N* = 1 in TV = 2. The standard scalar product in 
V is 

(u, v)v = / u ■ v + curl u • curl v + div u div v. 

Jq 

On the other hand, the natural space for (2) is 

JTi(n) = {u E H'iQ)” I u x n|r - 0} 
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with the scalar product 


(u, v)j = / u * v -f Vu • Vv. 

Jn 

We will show that while problem (1) has a unique solution in V, it may not 
be solvable in i/o t (f£), an d while problem (2) has a unique solution in 
the solution may not be unique in V . We will also discuss in what sense the 
divergence boundary condition in problem (2) is satisfied. 

It is shown in section 6 (under the additional assumptions that T is piecewise 
C 1,1 and the jump condition given in (6) is satisfied) that problem (2) (in 
#o t (f2)) is equivalent to: find u £ Hq^CI) such that 


1 

D> 

c 

II 

•-ta 

in ft, 

du 


— — • n + ku • n = 0 
on 

on r, 

u x n — 0 

on r. 


( 3 ) 


where k is the curvature of the boundary (k = divn), which is defined almost 
everywhere on T. This follows from the result (see Lemma 10 for the precise 
formulation) 

... du 

div u\ r = — • n| r + «u • n|r- 

It is convenient to begin with the weak formulations of problems (1) and (3). 
The next three sections deal with the coercivity of the corresponding weak forms. 
Following this we show the equivalence of the strong and weak problems and 
then examine the relationship of problems (1) and (3), and show the equivalence 
of problems (2) and (3). 


3 Weak Formulation in V 

Let Q and f be as above. To derive the weak formulation of (1) in V we substitute 
curl curl u for -Au in (1), multiply the first equation by v (belonging to K), 
the second by div v and integrate by parts using the boundary condition. For 
u, v £ V define 

a(u, v)= / curl u • curl v 4- div u div v. (4) 

Jn 

Then the statement of the weak problem is: find u £ V such that 

a(u,v)= / f • v Vv £ K (5) 

Jn 

The only hypothesis of the Lax-Milgram lemma which is nontrivial to check 
is the coercivity of a on F. It follows from the next compact embedding theorem. 
The embedding is actually a corollary of the regularity result in [1], where it is 
proved that V C We have provided a new and concise proof. 
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Theorem 1 V is compactly embedded in L 2 {{ T) N . 

The proof makes use of Murat’s div-curl lemma (see [8]). 

Lemma 1 (div-curl lemma). Let U be an open subset of H N . Let v n and 
w n be weakly convergent sequences in L 2 {U) N with the limits v and w corre- 
spondingly. Assume that ^divvn}^! and {curlWnJ^Lj lie in compact subsets 
of H~ l (U) and H~ l (U) N respectively. Then for every 0 e T)(U), 

I 0v n * w n — ► / 0 v - w as ti — ► oo. 

Ju Ju 

Proof of Theorem 1. Let u n be a weakly convergent sequence in V with the 
limit u. We want to show that u n converges strongly in L 2 (ft). The idea is 
to extend u n outside ft in two ways: v n and w n will be the extensions with 
“good” divergence and curl respectively. We construct the extensions so that 
supp v n Pi supp w n C ft and then use div-curl lemma for these sequences. 

Let U be an open ball containing ft. We start with the construction of v n 
(the extension with divv n in a compact set of H~ 1 (U)). For each n € N we 
define g n £ H l (U \ ft) by 


{ A g n = 0 in U \ ft, 
^ = u„ ■ n on dfl, 
g n — 0 on dU. 


Now define v n by 


and w n by 


Then 


and 


V n = 


u n in ft 
V 5 in in U \ O 


w 


div v, 


curl w. 


{ u n in ft 
0 in {7 \ ft. 

{ div u n in ft, 

0 in U \ ft 

{ curlu n in ft, 

0 in U \ ft 


so {divvn}^! and {curlw n }^ =1 are bounded subsets of L 2 (U) and L 2 (ft) N \ 
respectively (hence he in compact sets of H~ l (U) and H~ 1 (Q) N ). Moreover, 
v n — ^ v in L 2 (U) n and w n w in L 2 (U ) N , where v and w are similar 
extensions of u. 
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Now choose 0 G X>(17) such that 0 = 1 in ft. Using the div-curl lemma we 
get (note that v n ■ w n = 0 in U \Q) 

/ U 2 = / (j) V n • W n ► / <j)V • w — / u 2 . 

Jn Jc; it/ Jn 

This together with u n u in L 2 (S1) shows that u n — + u in L 2 (ft). 

We will use Peetre’s lemma (see [2], for example) in the following form 

Lemma 2 Let E, E\ and E 2 be Banach spaces . Let A\ and A 2 be continuous 
linear operators from E to E\ and from E to E 2 respectively. Assume that there 
exists C > 0 such that 

IMIe < C(\\Aiu\\e 1 + \\A2u\\e 2 ) Vu G E. 

Assume also that KevAi = {0} and that A 2 is compact. Then there is C\ > 0 
such that 

IMIe < Ci||-4iix||£;i Vu g 
T he coercivity result is then 
Lemma 3 There is C > 0, independent of v, such that 

a(v, v) > C||v||y Vv G V 

Proof We use Lemma 2 with E = V, E\ = L 2 (Q) x L 2 (ffc) N , E 2 = L 2 (Q) N 
and ^4iu = (div u, curlu), A 2 M — u. By Lemma 1 the operator A 2 is compact. 
We must show that Ker A\ = {0}. But 

u G Ker A\ <=> div u = 0 and curl u — 0 in $7, and u x n - 0 on F. 

This implies that u = 0 (by using the gradient potential, for example). So by 
Lemma 2 the result follows. 

The Lax-Milgram lemma now implies that (5) has a unique, stable solution 
in V . 

4 Weak Formulation in H^Q) 

For deriving the weak form of (3) we shall make additional assumptions about 
the smoothness of the boundary. In the following Q is a bounded subset of , 
N = 2 or 3, with a connected, Lipchitz and piecewise C 1 ’ 1 boundary T, i.e. 

T = U j =1 Tj with Tj H Tt= 0 for i ^ j, and Tj, j — 1 , . . . , n are of class C 1,1 . 
Let 7 denote the set of “edges and corners” , i.e. 

n 

7= U (f«nr,). 
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We assume that the jump of the normal on 7 is bounded below, i.e. there is 
6 > 0 such that for all x E 7 with x E f { n f j , i ^ j we have 

n^x) • n J (x) <1 — 6 , (6) 

where n* (x) and n J (x) are the limits of the unit outer normals when approaching 
x from and Tj, correspondingly. 

This condition is satisfied for all polygons and polyhedra. Also included are 
polygons and polyhedra with curved sides. It excludes three-dimensional bodies 
which have points like the tip of a cone, and bodies with edges which “flatten 
out”, i.e. the angle between the faces gets arbitrarily close to 7r. This condition 
is needed for proving some results about trace operators in the next section, 
and implicitly in the proof of Lemma 4, which gives the relation between the 
two weak forms. 

To obtain a suitable weak form for (3) we proceed formally, multiplying the 
equation by a test function v e and integrating by parts as usual. We 

get 

[ Vu • Vv - f ^ • v = f f v. 

Jn Jr dn J Q 

Rewriting the boundary term as a sum over the smooth boundary pieces and 
using the boundary conditions on u and v gives the weak problem: find u € 
Hg t (Q) such that 

[ Vu-Vv + V f /cu ■ v = [ f • v VvEi/o\(^)- (7) 

Jn j =1 Jrj Jrt 

Since k E L°°(rj), by the usual trace theorem the boundary term is well de- 
fined. Note that, in this framework, the divergence boundary condition (or the 
equivalent one in problem (3)) is natural. 

Proving coercivity of the weak form (7) is a nontrivial matter because the 
curvature can be of either sign. The result will follow from the next lemma 
whose proof is deferred to the following section. 


Lemma 4 Assume that satisfies the assumptions made in the beginning of 
the section. For u, v E we have 


where Tj, 
in (4). 


a(u, v) = / Vu ■ Vv + / ku • v, 


(8) 


j = 1,.. . ,n are the smooth (C 1 ' 1 ) pieces ofT, and a is the same as 


This means that the weak problem (7) is in fact equivalent to: Find u E 
HQ t (Q) such that 


a(u, 



f v 


Vv € 


( 9 ) 
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in which the equation is the same as in (5) but the spaces are different. 

In fact, one can get (9) from problem (1), but they are not equivalent: 
in general we cannot assume that the solution of (1) belongs to i/o t (S7), and 
we cannot prove later that the solution of (9) satisfies divu = 0. We can 
also formally get (9) from (2) by substituting — Au by curl curl u — Vdivu, 
multiplying by v and integrating by parts using the boundary conditions, but 
this is not easy to justify rigorously (for the solution of (2), curl curl u and 
Vdivu may not be in L 2 ). 

The coercivity result for (9) is: 

Lemma 5 There is C > 0 such that 

°( u ) v ) > C||u|| ! Vu€f&(n). 

Proof. We use Lemma 2 with E = //^(fi), E\ = L 2 (fl) x L 2 (Q) N , E 2 = 
L 2 (Fl) N x L 2 (F) n , Aiu = (divu, curl u) and A 2 U = (u, u|p). The operator 
A 2 : E — > E 2 is compact because it is bounded as an operator from E to 
/f 1 (Q) x H l / 2 (T) N and the latter is compactly embedded into E 2 . We already 
showed (in the proof of Lemma 3) that Ker A\ = {0}. The result follows. 

Applying the Lax-Milgram lemma to the weak form (9) (or, equivalently, to 
(9)) and using the coercivity result above shows the existence of a unique stable 
solution to (9). 

Note that to prove the coercivity of a in H$ t (f£) we do not need the additional 
smoothness assumptions made in the beginning of this section: they are needed 
to make sense of the term containing the curvature in (7) and to show later that 
the divergence boundary condition is satisfied for the solution of problem (9). 
They are not needed for uniqueness of solution of (9). 


5 Proofs for previous section 


This section will present the proof of Lemma 4. We will assume throughout the 
section that the assumptions made in the beginning of the previous section are 
satisfied. 

First we need some preliminary results. We begin by recalling some facts 
about the trace spaces. First, for Lipchitz T C R N , the following is an equivalent 
norm on H l f 2 {F) (see e.g. [3]): 

m 2 HU Hr) = ||*|ll.(r) + H 


We will also need the spaces Hqq 2 (Fj), j = 1, . . . ,n, which consist of all func- 
tions in H l / 2 (Tj) whose extension by zero to F belongs to H l / 2 (F). The norm 
of a function in Hq q 2 {Fj) is the norm of its extension by zero in H l / 2 (F), i.e. 


,2 11 j 1 , 2 , f f l<£(x)- 0 (y )| 2 JO 

+ JrJr, |x-y|" ’ y 

+2 [ f 

J r\r 3 Jr D l x _ y 


w dS-xdSy. 
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1/2 

We will use the same notation for a function in H 0 q (1^) and its extension by 
zero to r. 

We will also need the following result: 


Lemma 6 The trace operator v i— ► v ■ is a continuous linear operator from 
Him onto tfoo 2 ( r i). 


Proof. We use the fact that the usual trace operator v m v|r is a continuous 
linear operator from onto i/ 1 - / 2 (r) Ar . First we will show that if v x n = 0 

and the jump condition (6) is satisfied, then ||v • n|| H i/ 2 ^ r ^ < C||v|| H i/ 2 ( r ). Let 

cj) = v • n on Tj and <p = 0 on F \ Tj. Clearly \\<f>\\H 1 / 2 (r j ) < ||v||jyi/ 2 (p >, so we 
have to show only that 

L, L 

Because of the jump condition we can find e > 0 such that 

X € r,-, y £ r \ Tj, |x — y| < e =s> n(x) ■ n(y) < 1 - 

For x 6 Tj y y G T \ Tj such that |x — y| < e we have 
|v(x) - v(y)| 2 = |v(x)| 2 + |v(y)| 2 - 2(v(x) ■ n(x))(v(y) • n(y))(n(x) ■ n(y)) 

> |0(x)| 2 + |v(y)| 2 - 2|4(x)| |v(y)| (l - 0 

> I0(x)| 2 ^1 ~ ^ = M x )| 2 (<* - j) ■ 


5 2 

Since S — — > 0 (note that 0 < S < 2, otherwise the jump condition could not 
be satisfied), we have 

W x )| 2 


/ / 

Jr\r, Jr, l x - , 

-/ / 

Jr\r j J r d \ 


ijy dS x dSy 


W x )| 2 


dS x dS v -j- 


7jy u ^x u ^y 


Ir\ r, Ir , 


l<A(x )| 2 


r\rj Jr,nB'( y ) 


i\s.( y> l x -yl ; 

<4 / / |<£(x)| 2 dS x dS y + -V / / 

€ Jt\t 3 Jr, o - T Jr\r, Jr, 

— ^ll v l!H 1 / 3 (r)* 


|x-y| 


jy dSxdSy 


l v (x) — v(y)| 2 

| x _ y | N 


Consequently, the operator v i— > v ■ n|rv is a continuous linear operator from 
Him to tfoo / 2 (r,). 

Let <f> G be given. Extend it by 0 to T. Then <j>n € H l / 2 (T) and 

we can find v e such that v|p = <f>n and ||v||i < C\\(f>n\\ H i/ 2 ^ < 
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Ci|| 0 || „i/2> r v Clearly v £ HQ t (tt) y so the trace operator is onto. 

“oo ^ 3 ) 

The next lemma gives a formula for calculating the divergence on the bound- 
ary. 


Lemma 7 For smooth u the restriction of div u to Tj satisfies 


r I /j- M 

div u r = I divpu T + acu • n + — — - n 

V an ; | r . 

where u T is the tangential component of u and divp is the divergence of a tan- 
gential vector field in the tangential coordinate system. 

Proof. To calculate div u at some xo 6 Tj let us fix a coordinate system i, j, 
k (if N = 2 then omit j), where i, j are tangent to T at xo and k is the normal at 
xo (we can do this for every xo in the interior of Tj). Denote u r = u — (u ■ n)n. 
Then 

i / \ <9u . . du du 

divu(xo) = -^r(xo) * i+ ^j-(xo) *J + ^( x o) * k 


du r 


du r 




a[(u ■ n)n] a[(u ■ n)n 


di 


(xo) ■ i + 


dj 


-(xo) - j 


du. . . 

+ g-(xo) • n(x 0 ) 


a u 


(div r u T + /cun+^ 


") 


where 


div r u T (x 0 ) = ^(xo) • i + ^ L ( x o) ' J- 


Now we can prove lemma 4: 

Proof of Lemma 4 • For v € #Q t (f2) and smooth u (we do not require 
u x n = 0 yet, because it does not seem straightforward to show the density of 
smooth functions with zero tangential trace in HQ t (fl)) we have 

/ curl u • curl v + div u div v = / curl curl u • v — Vdivu • v + / (divu)v n 

Jn Jn Jr 

= — / Au • v + / (div u)v • n 

Jn Jr 

— J Vu ■ Vv + J ^divu — • nj v • n. 


By Lemma 7 


/ curl u • curl v + div u div v = / Vu ■ Vv + V' / (divpu r -b «u ■ n) v ■ n. 

Jn Jn j= l Jr i (10) 
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For any 0 G HQQ 2 (Fj) (extended by zero to F) we can find v € suc ^ 

that v|r = <j>n and ||v||i < C||0||„i/2, r v From (10) we get 

/ (div r u T + KU • n) <f> < C'llvllillulli < C"||0||^/2 (rj) ||u||i 

JTj 00 

V0 € tfoo 2 ( r j)> Vu smooth. 

By density of smooth functions in f/ 1 (J7) and Hahn-Banach theorem we can 
continue the mapping u ♦— > (divru r + ku ■ n)|r . uniquely to a continuous linear 

mapping from H l (Fl) N to (HQQ 2 (Fj)y. Since for any v G HQ t (Fl), v • n G 
H^Tj) (by Lemma 6), (10) holds for all u € H 1 ^)" and v € (the 

integrals over Fj have a meaning as a duality pairing between (HfJ 2 (Y j))' and 
^oo 2 ( r l))- 

To finish the proof we need to show only that if u G i/o t (f2), then divpu r = 0 
on Fj. For this first note that since uh ku • n\r j is linear and continuous from 

H 1 (Fl) N to L 2 (Fj) C ( HQQ 2 (Fj )) f , the operator u i— ► divru r |iv is also continuous 
from H l (Q) N to (Hqq 2 (F j))' , and then use the fact that 

<div r u T ,0) rj = ~ f u T * V r </> V<j> G Cl(Fj) 
and the density of C \ (r j ) in Hqq 2 {F j). 

6 Equivalence of the Strong and Weak Forms 

In this section we will prove that (1) in V is equivalent to the weak problem 
(5), and that (3) in i/o t (fl) is equivalent to the weak problem (7). It follows 
that problem (1) is uniquely solvable in V and problem (3) is uniquely solvable 
in i/o t (H). We will also show that problems (2) and (3) are equivalent with the 
divergence boundary condition having a meaning in the sense of traces. 

Let us first deal with problem (1). 

Lemma 8 Let and f be as m Section 2. Then the problem of finding u G V 
satisfying (1) is equivalent to solving the weak problem (5). 

Proof. Let u G V be a solution of (1). Then — Au — curl curl u G L 2 (ft)) 
and so the formal calculations leading to (5) can be rigorously justified. 

Conversely, if u G V satisfies (5), then by using v G P(S7) N and transferring 
all derivatives to v we get — Au = f in the sense of distributions. For any 
<j> G 'D(Q) we can find g G Hq(FI) such that A g = f>. Now use v = V# in the 
weak form (note that V<? G V). We get 

[ (div u )<f> = [ i Vg = - ( (div f )g 4- [ (f • n )g = 0, 

Jn Jn Jn Jr 
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so divu = 0 in the sense of distributions. Consequently the problem (1) (in V) 
and the weak form (5) are equivalent. 

Corollary 1 The problem (1) has a unique solution in V . 

In the proof of the theorem we used the fact that Vg £ V. If we were dealing 
with the weak form in then Vg may not have been a legitimate test 

function (in the case g £ H 2 (i 7)), which would make it impossible to prove that 
divu = 0. This is the difference between the weak problems (5) and (9), even 
though they appear very similar. 

The corresponding proof for the problem (3) is a little more difficult, but 
follows the same pattern. 

Lemma 9 Let ft and f be as in Section 2. Assume in addition that the boundary 
T is piecewise C 1,1 and the jump condition (6) is satisfied. Then the problem of 
finding u £ 77^(17) satisfying (3) is equivalent to solving the weak problem (7). 

Proof We first have to show that the boundary condition 

du nr 

— • n -f ku • n = 0 on 1 
an 

has a meaning for u £ 77^ ($7) with Au £ L 2 (S7). Since Vu £ L 2 (f2) iV and 

du 

divVu £ Ir(fJ), the gradient has a normal trace on the boundary — £ 

H~ 1/2 {F) (see [2] for example). For <j> £ H^ 2 (Tj) we have <f > n £ H l ^ 2 {F) 
(using the extension of (p by zero) and therefore the boundary condition may be 
interpreted as 

• + ku, </>n) r = 0 V^fl Hli 2 {T 3 ), (11) 

j= 1 


or, equivalently, 

+ «u, v) r = o Vv € H&n), (12) 

where (•, -)r is the duality pairing between and i/ 1 / 2 (T). 

Let u £ H^ t {Q) be a solution of (3). Using (12) we can justify to formal cal- 
culations leading to the weak form (7) (the boundary integrals in the weak form 
should be understood as duality pairings between (Hqq 2 (F j)) f and HoQ 2 {Fj)). 

Conversely, if u £ satisfies (7), then by using v £ V(Fl) N we get as 

before — Au = f in the sense of distributions. For any v £ #^(£1), integrating 
by parts in the weak form (note that Au £ L 2 (£l) N ) we get 


-/ 

Jo 


Au • v + 


(- 

'an 


,v)r+yf 




KU • V 


In 


f v, 


hence (12) is satisfied. Consequently the problem (3) (in 7/^(0)) and the weak 
form (7) are equivalent. 
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Corollary 2 The problem (3) has a unique solution u in HQ t (fl). 

Now we will show that problems (2) and (3) are equivalent. The only dif- 
ference in these two problems is in the boundary condition, so we have to show 
that the divergence boundary condition makes sense and is equivalent to the 
boundary conditon in (3). This is done in the following lemma. 

Lemma 10 Assume f] satisfies the conditions of Lemma 9. For u £ HQ t (fl) 
with Au € L 2 (ft) N , and any j £ {1, . . . ,n), the trace of divu on Tj exists and 
belongs to (H^ 2 (T j))' , and 

,, du 

divu|rv = ^ ■ n| ri + «u ■ n|ry 

Proof This follows from Lemma (7) which gives the formula for smooth u 
without the condition uxn = 0. By density this can be extended to u £ H 1 (f 1) N 
with Au £ L 2 (fl) N . Then use the fact that for u £ .ffo t (S7) we have div pu T = 0 
(see the end of proof of Lemma 4). 

Corollary 3 Problems (2) and (3) are equivalent. 

Let us summarize the results in a theorem: 

Theorem 2 Assume that ft is a bounded subset of H N with a connected, Lip- 
schitz boundary. Let f £ L 2 (ft) N with divf = 0 be given. Then 

a) The problem (1) in V is equivalent to the weak form (5) and has a unique 
solution uy £ V . 

b) If in addition the boundary is piecewise C 1,1 and satisfies the jump con- 
dition (6), then the problems (2) and (3) in HQ t (fl) are both equivalent to weak 
problems (7) and (9) y and have a unique solution u h G Ho 

7 Relationship of the V and HQ t (Q) Formulations 

Let uy £ V be the solution of (5) (or the strong form (1 with the interior 
divergence condition) and let u h € #Q t (f2) be the solution of (9) (which under 
the additional smoothness conditions is equivalent to problems (2) and (3) with 
the divergence boundary condition). In this section we will discuss when the 
two solutions are the same. The following is clearly true. 

Lemma 11 Let and f be as in Section 2. Assume in addition that the jump 
condition (6) is satisfied. Then the following are equivalent 

a) uy = u h; 

b) divu h —0 in fl; 

c) Uv € 

Let us now examine when uy — u h. First, note that lemma 5 implies 
that HQ t (fl) is closed in V. We have two possibilities: either V = HQ t (fl) 
or F / HQ t (fl). In the first case the solutions are obviously the same (the 
corresponding weak forms are exactly the same). The following lemma gives 
the necessary and sufficient conditions for the equality of the spaces. 
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Lemma 12 Assume that SI is a bounded subset of R N with a connected Lip- 
schitz boundary. The spaces V and HQ t (Sl) are equal if and only if the homo- 
geneous Dirichlet problem for the scalar Poisson equation with any right hand 
side in T>( SI) has a solution in H 2 (Sl). 

Proof. Suppose that the regularity condition holds. For any given w £ V 
we can solve the weak problem: find u £ HQ t (Sl) such that 

a(u, v) = a(w, v) Vv £ i/J f (ft) (13) 

(i.e. u is a projection of w onto HQ t (Sl) corresponding to the scalar product 
a(-, •) ). As in the proof of Lemma 8 we can show that A(u — w) = 0, and 
using the //^-regularity of the solution of the Dirichlet problem for the Poisson 
equation we also get div (u — w) = 0. Since (u — w) x n|r = 0, by the uniqueness 
of solution of (1) we have u = w, i.e. w £ HQ t (Sl). Consequently V = Hq ( (Si). 

Now suppose that there is g £ Hq (SI) such that A g £ T>(Sl), but g £ H 2 (Sl). 
Then S/geV, hut so V^H^ t (Sl). 

Note that the lemma remains true if one substitutes L 2 (Sl) for V(St). 

Corollary 4 If SI is convex or has a C 1 ' 1 boundary then V = HQ t (Sl). 

This follows from the classical regularity theory about the smoothness of the 
solution of the Dirichlet problem for the Poisson equation. It is actually enough 
to assume that Si is locally convex near the points where the boundary is not 
C 1,1 (i.e. these points have a neighborhood whose intersection with SI is convex). 

The question of the equality of these spaces (and analogous ones with a 
boundary condition for the normal trace) has been partly answered by different 
authors, but we are unaware of results as sharp as these presented here. In the 
classical book [2] it is proved that V = HQ t (Sl) if SI is a convex polygon or has 
a C 1,1 boundary. Necessary and sufficient conditions for V = HQ t (Sl) to hold in 
R 2 are given in [6], namely that the scalar Poisson equation with L 2 (Sl) right 
hand side must have a solution in H 2 (Sl) both with Dirichlet and Neumann 
condition. As we saw, the regularity of solution of the Neumann problem is 
not needed, and in R 2 it actually follows from the regularity of solution of the 
Dirichlet problem. 

In the other case, when //^(O) and V are not identical, since HQ t (Sl) is 
closed in V, there is f £ L 2 (Sl) with divf = 0 for which the solution of (1) is in 
V, but not in HQ t (Sl). This means that the interior and boundary divergence 
formulations are not equivalent and give different solutions for this f . We will 
show how to construct f for which u y ^ u h whenever V ^ HQ t (Sl). 

Lemma 13 IfV HQ t (Sl) than there is f £ L 2 (Sl) with divf = 0 such that the 
corresponding solution of (1) is not in HQ t (Sl). 

Proof If V ^ i^oi(^) then we can choose a nonzero w £ V such that 

a(w,v) = 0 Vve (14) 
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First note that curl w ^ 0: indeed, if we had curlw — 0 then 

f divw-divv = 0 Vv € #o t (R), (15) 

Ju 

and hence (by using v € £>(0)^)), Vdivw = 0, or divw =const. It follows 
from (15) that div v = 0 for all v € //^(fi), which is a contradiction. 

Now let f £ V be a solution of 


a(f 


,v) = / 

Jn 


curl w • curl v Vv € V. 


(16) 


Note that f ^ 0 and divf = 0. We claim that for this f problem (1) has a 
solution uy E V which is not in Indeed, using the weak problem (5) 

and (16) we get 



curl u • curl w = 


«( u,f)= f If I 2 ^ o, 

Jn 


which by (14) implies that u y ^ HQ t (fl). 


We can also give a more tangible example of a solution of (1) which does 
not belong to //^(Q). Suppose $1 C R 3 is such that in a neighborhood of some 
x 0 G T its boundary consists of 2 planes with the interior angle a > tt. Choose 
a cylindrical coordinate system (r, 0 , 2 ) with the origin at xo, 2 -axis along the 
edge and 0 = 0, 0 = aon the boundary. Let 0 £ C°°(f2), depending on r and 
2 only, be such that 0 = 1 near the origin and 0 = 0 outside a neighborhood of 
the origin where the pieces of T are planar. Put 

u(x, z) = curl ^0(r, z)r n ^ a cos — , 

where e z is the unit vector in direction of 2 -axis. Then div u = 0 and 

uxn=H 0(r, z)r n ^ a cos — J e z = 0 on I\ 
on \ a J 

It follows that u is a solution of (1) with f = -Au E L 2 (Q) and divf = 0. In 
addition, f = 0 near the origin. Since u £ H^ t (Q) and the problem (2) with the 
same f has a solution u £ i/o t (fl), we must have u ^ u and divu ^ 0. 

This example is essentially two-dimensional. To make it work for Q C R 2 
in a neighborhood of a corner with the interior angle a > 7T we just use 0 
independent of z. Note that the third component of u is zero, and we get 
f £ C°°(Q) 2 . 

We will state these results as a theorem: 

Theorem 3 Assume that Q is a bounded subset of H N with a Lipschitz bound- 
ary. The following are equivalent: 

a) uy = u h for any f € L 2 (fl) N with divf = 0; 

b) V= 

c) the Dirichlet problem for the scalar Poisson equation with any smooth 
right hand side has a solution m H 2 (Sl). 
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8 Generalization to Helmholtz equation 

Most of the arguments above remain true in case of Helmholtz equation. The 
only difference is in the coercivity, which now will not be sufficient for using 
Lax-Milgram. Let us sketch the results briefly. The analog of problem (1) is 
now: find u G V such that 

— Au — k 2 u — f in 

i div u = 0 in Q, (17) 

u x n = 0 on T. 

We assume that f G L 2 (fl) N and the compatibility condition divf = 0 is satis- 
fied. To formulate the weak problem let us define 

b( u,v) = / curl u ■ curl v + div u div v — fc 2 u • v Vu,vgV. 

Jn 

Then the weak problem corresponding to (17) is to find u G V such that 

b( u,v) = f f v Vv € V. (18) 

Jn 

The coercivity result is that b(u, u) > C||u||y — A; 2 ||u|| 2 2 (if k 2 is complex, 
we should use the real part). We cannot use Lax-Milgram directly. Instead 
we can use the Fredholm alternative to show the existence and uniqueness of 
the solution. Define operator Ay : L 2 ( f!) — > L 2 (Q) by the requirement that 
Ayg G V is the solution of a( u, v ) = / fi g*vVv6V r (note that we do not need 
the compatibility condition div g = 0 to guarantee the solvability of the weak 
problem). By the Lax-Milgram Lemma and the compact embedding theorem 
(Theorem 1) Ay is a compact operator. Then the solution of (18) is the solution 
of u = Ay{k 2 u 4- f), and if 1/A: 2 is not an eigenvalue of Ay then (18) has a 
unique solution u G V. The operator Ay has two kinds of eigenfunctions: ones 
that are divergence-free, and others which are gradients of the eigenfunctions of 
the scalar laplacian with Dirichlet boundary conditions. Assuming that 1/k 2 is 
not an eigenvalue of Ay , one can prove the equivalence of the strong and weak 
problems similarly to the proof of Lemma 8 (in proving that the solution of the 
weak form satisfies the divergence condition we need solvability of the scalar 
Helmholtz equation with Dirichlet boundary conditions, which follows from the 
fact that if 1/A; 2 is not an eigenvalue of Ay, then it is not an eigenvalue of the 
scalar laplacian with the Dirichlet boundary conditions) . 

The analog of problem (2) is to find u G Hq^SI) such that 

r — Au — A: 2 u = f in O, 

< div u = 0 on F, (19) 

u x n = 0 on T 

V 

and the corresponding weak problem is to find u G HQ t (£l) such that 

f >(u,v)= f f V Vvg4(9). (20) 

Jn 
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One can deal with the coercivity as above, using the operator Ah : L 2 ( f2) — ► 
L 2 (S7) defined by: Ang 6 i7o t (f2) is the solution of a (u, v) = Jg-vVv € i/o t (D). 
Equivalence of the strong and weak forms follows similarly to the zero frequency 
case. 

The question when the two weak problems (18) and (20) give the same 
solution can be answered as in section 7. Most of the results presented there do 
not depend on the frequency. The counterexamples are also easy to modify: if 
u v £ HQ t (£l) is a solution of (1), then it is also a solution of (17) with the right 
hand side f — k 2 uy. Theorem 3 remains true if we assume in addition that 1/A: 2 
is not an eigenvalue of either Ay or Ah and let u y and u h be the solutions of 
the weak problems (18) and (20) correspondingly. 

Similar results can be proved for time-dependent problems, such as vec- 
tor wave equation and Maxwell’s equations. Even with divergence- free initial 
conditions, using a weak form in 77^ (ft) (similar to (20)) may give a solution 
not satisfying the divergence constraint. Again, this can happen only if the 
i7 2 -regularity of the Dirichlet problem for the scalar Poisson equation fails. 


9 Conclusions 

We saw that in some cases (e.g. for nonconvex polygons) the boundary diver- 
gence formulation is not equivalent to the interior divergence formulation, and 
the weak form (9) or (20) (the weak forms in 77o*(ft)), while uniquely solvable, 
may not have a divergence-free solution. In this case all approximation meth- 
ods based on these weak formulations (e.g. the usual finite element method, 
the least squares method) will converge to a spurious solution (not satisfying 
the divergence constraint). To avoid the spurious solution, one should use the 
weak form in V , but this is not straightforward: one has to use basis functions 
which are in V , but not in ifQ t (ft) (e.g. the singular solutions around the re- 
entrant corners), in addition to the usual ones. The same is true about the 
penalty method, where the weak form contains a penalty parameter in front 
of the term with divergence. If the solution of the original problem is not in 
77o t (ft), then any choice of the penalty parameter will result in a spurious solu- 
tion — in fact, the smaller one makes the divergence (by choosing larger values 
of the parameter), the bigger the error is in the curl of the solution. 

On the other hand, when the formulations are equivalent (this is the case 
when the solution of the Dirichlet problem for the scalar Poisson equation is in 
H 2 ( ft)), then one can use the weak forms in i/Q t (ft). This enables one to use 
simpler algorithms, e.g. finite element methods with piecewise linear test and 
trial functions. The equivalent formulation in terms of gradients (as in (7)) may 
be especially useful, since it decouples the field components inside the domain. 
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